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BCS Diquark Condensation in the 3 + Id Lattice NJL Model 

David N. Walters 

'^Department of Physics, University of Wales Swansea, 
Singleton Park, Swansea SA2 8PP, United Kingdom. 

We present preliminary evidence of BCS diquark condensation in the 3+1 dimensional Nambu-Jona-Lasinio 
(NJL) model at non-zero chemical potential {jj.) on the lattice. Large results are used to match the model's 
parameters to low energy, zero density phenomenology. A diquark source j is added in a partially quenched 
approximation to enable the measurement of lattice diquark observables. In particular measurements are made 
of the diquark condensate and susceptibilities as functions of j which support the existence of a BCS phase at 
high /i. 
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1. Introduction 

Despite recent advances in the study of QCD 
at small but non-zero density on the lattice 0| the 
sign problem still persists in the study of dense, 
low temperature, exotic matter phases. One ap- 
proach is to study four-fermi theories which, un- 
like SU(2)c' QCD (in which the qq baryons are 
bosonic), can display Fcrmi-surface effects such as 
BCS superfluidity. Although these models exhibit 
no confinement, one variant, the SU(2)l x SU(2)fl 
NJL model, possesses the same global symmetries 
as 2 flavour QCD, which should allow us to inves- 
tigate the phase structure of this theory. 

2. Lattice Model and Parameter Choice 

The model studied here, which uses staggered 
lattice fermions, is the 3+1 dimensional version 
of that studied in Q and references therein. In 
particular the model has the action 
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where the bispinor is written in terms of in- 
dependent isospinor fermionic fields via — 
(x, X*'^); s-nd the auxiliary bosonic fields a and 
TT are introduced in the standard way. Written in 
the Gor'kov basis the fermion matrix is 
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where the matrix M is identical to the 3 + Id form 
of that given in H. The diquark sources j and 
J differ from those in by a factor of 2 which 
allows us to identify them with Majorana masses. 

Being a theory of point-like interactions the 
NJL model has no continuum limit for c? > 4 
leaving the physics sensitive to simulation param- 
eters. In order to ensure we simulate in a physical 
regime these parameters are fitted to physical ob- 
servables calculated in the large N limit at /z = 
1^ . In particular we evaluate dimensionless ratios 
between the pion decay rate /tt, the constituent 
quark mass m* and the pion mass mir, and fit to 
the phcnomcnological values of 93MeV, 350MeV 
and 138MeV respectively. Using these fits we de- 
termine the bare quark mass amo = 0.002 and 
the inverse coupling a^/g"^ = 0.565. We also ex- 
tract the large N lattice spacing a^^ = 1076MeV, 
which allows us to present results in physical 
units. 



3. Preliminary Results 

The above model was simulated with a Hy- 
brid Monte Carlo algorithm on x Lt lattices 
with Ls = Lt = 12, 16 and 20. A partially 
quenched approximation was used in which the 
diquark sources are set to zero during the update 
of the auxiliary fields but are non-zero during the 
measurement of diquark observables. 

In order to study chiral symmetry breaking we 



2 



measure the chiral condensate (xx) ^^'^ baryon 
number density ns, which are defined by 
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Figure 1. Chiral condensate (xx) ^ind baryon 
number density as functions of /it with lattice 
results compared to the large N limit. 



Plots of (xx) and as functions of chemi- 
cal potential are presented in Figure |l|. Both the 
large N solutions and the lattice data have been 
extrapolated to the limit V^^ 0. Chiral sym- 
metry is broken in the low fi phase, as signalled 
by an non-zero chiral condensate, and is approx- 
imately restored as ^ is increased through what 
appears to be a crossover. The order of the chi- 
ral phase transition in the NJL model is strongly 
dependent on the parameters used [1). 

In the diquark sector, the operators 

qq±{x) = x**^ ±x jx*'^(a;) (5) 

allow one to define the diquark condensate as 
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Figure 2. Diquark condensate vs. j for high and 
low fi. 



where j± = j ± j. 

Figure || shows the condensate plotted as a 
function of j extrapolated to L^^ — > (we sim- 
ulate with j = '] throughout). For fi — Q a, 
quadratic fit through the data is consistent with 
the condensate vanishing as j 0. For high /i the 
fit is consistent with a non-zero diquark conden- 
sate which is similar in magnitude to the chiral 
condensate in the vacuum. We believe that finite 
volume effects cause the points with j < 0.2 to 
fall below the curve. In particular, in the high 
fi phase as J — > we expect an exact Goldstone 
mode which leads the correlation length to di- 
verge. For this reason these points have been ig- 
nored during the fits. 

Finally, if we define susceptibilities 

X± =E(99±(0)'?'?±(^))' (7) 

X 

we can clearly distinguish between the two phases 
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by studying the ratio 



(8) 



If a U(l) baryon number symmetry is manifest, 
the two susceptibihties should be identical up to 
a sign factor and the ratio should equal 1. If the 
symmetry is broken, the Ward identity 



X-\j_=o 



(gg+) 
j+ 



(9) 



predicts that x_ should diverge and the ratio 
should vanish. 
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Figure 3. Susceptibility ratio R vs. j for high and 
low fi. 



Figure || shows the susceptibility ratio as a 
function of j, again extrapolated to L^^ 0. 
A linear fit through the /i = data is consistent 



with a manifest U(1)b symmetry. For fia = 1.0 
we see markedly different behaviour, with a fit 
(again discarding j < 0.2) suggesting an inter- 
cept close to zero and therefore a broken baryon 
number symmetry. 

4. Summary and Outlook 

We have provided preliminary evidence of a 
BCS superfluid phase in the 3-|-ld NJL model at 
high /i, with an order parameter approximately 
equal in magnitude to the chiral order parameter 
aX II — 0. The behaviour of the susceptibility ra- 
tio defined in further supports a broken U(1)b 
symmetry in the high fi phase analogous to the 
broken \J{1) em in BCS superconductors. This 
study provides the first evidence of a BCS phase 
in lattice field theory, which we interpret as evi- 
dence of a colour-superconducting phase at high 
fj. in QCD. 

Clearly, more data are required for intermedi- 
ate values of /i to understand the nature of the 
diquark transition. To have more control over 
the finite volume effects shown in Figures |^ and 
H it may also be necessary to study lattices with 
Ls ^ Lt- Finally, in the future we wish to study 
the fermion dispersion relation. This will allow us 
to present a comparison between the chiral mass 
gap S and the BCS superfluid gap A which can, 
in principle, be experimentally determined. 
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